BARI-TH 519/05 

Smeared gap equations in crystalline color superconductivity 

M. Ruggieri 

Universitd di Bari, 1-70126 Bari, Italia 
and 

I.N.F.N., Sezione di Bari, 
1-70126 Ban, Italia 

Abstract 

In the framework of HDET, we discuss an averaging procedure of the NJL quark-quark interac- 
tion lagrangian, treated in the mean field approximation, for the two flavor LOFF phase of QCD. 
This procedure gives results which are valid in domains where Ginzburg-Landau results may be 
questionable. We compute and compare the free energy for different LOFF crystalline structures. 
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I. INTRODUCTION AND THE LOFF STATE 



The behaviour of QCD at very high baryon density and low temperature has recently 
attracted a lot of interest. In these conditions, quarks are expected to deconfine [jj and 
occupy (in momentum space) large Fermi spheres. At very high densities the relevant inter- 
action is the one gluon exchange, which is attractive in the antisymmetric 3 color channel. 



Thus, a Cooper- 
tivity 



£e pairing phenomenon is expected to occur. This is color superconduc- 



a Uooper-nke pairing pnenomenon is expe< 
(see H 0, Q, Q 



14j for reviews). 



Regimes of high baryon densities and low temperatures are expected to be realized in the 
core of compact stars; therefore these superdense objects could be the places where color 
superconductivity is realized in nature. Apart the astrophysical speculations, the study 
of color-superconductive QCD is an intriguing challenge itself because its understanding 
implies a deeper understanding of the QCD phase diagram. 

In nature, as a consequence of electrical and color neutrality and of the different masses, 
the Fermi momenta of the quarks should depend on their color and their flavor. When the 
difference of the Fermi momenta of the pairing quarks are too different, a superconductive 
ground state in which the Cooper pairs have a net momentum is energetically favored with 
respect to the usual, zero momentum state. The resulting state is known as LOFF phase, 
and has been studied in the sixties in the context of condensed matter physics 



F phase, 

HQ. 



Possible realizations of the LOFF state in the frame of color superconductivity have been 
considered for the first time in Refs. G, 13, Q (see Q for a review and Refs.fl fl Q 
for recent developments). 

A color-superconductive phase is characterized by a nonzero expectation value of a bi- 
linear quark operator, namely a condensate. In the LOFF phase with two flavors (u and d) 
one has 

(^C 7 5^)ocA(r) e ^ 3 e Jj3 , (1) 
where a, (3 are color indices and i, j denote the flavor. A three flavor case (u, d and s) has 



been recently considered in Ref. 2J|. As for the space dependence of the condensate one 
usually decomposes A(r) as a sum of plane waves. The most simple ansatz is [lfij ]: 

A(r) = A e 2iqr , (2) 

and the corresponding superconductive phase is known as FF state. The net momentum 
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of the Cooper pair is 2q, which is the same for all the pairs while its direction is chosen 
spontaneously. 

One of the most important problem is to calculate the gap parameter A in Eq. (J2J). For 
the FF state this problem can be solved exactly: one has simply to shift the momenta of 
the paired quarks, p u — > p u + q, — > pd — q so the q dependence of the gap parameter 
disappears. Then one can write a self-consistency and exactly soluble equation for A. In the 
context of QCD this procedure has been applied in Ref. where a diagrammatic approach 
to the FF color superconductive state is considered. Unfortunately, this procedure cannot 
be applied to the general case of the linear combination of iV plane waves. To overcome 
this difficulty, one can employ a Ginzburg-Landau (GL) expansion of the general LOFF free 
energy functional Q. GL expansion works where A is small when compared to the typical 
mass scales of the model, in this case q and Sfi. In the framework of the GL expansion, in 
Ref. 13| it is conjectured that a face centered cube is a good candidate for the crystalline 
color superconductive state. 

In this paper we wish to define a procedure that allows to study the crystalline color 
superconductor with a generic crystalline structure in a space of parameters which is com- 
plementary to the GL one. This procedure is based on a weighted average of the interaction 
lagrangian over the lattice cell. The paper is organized as follows: in Sec. |H] we discuss 
the smearing procedure used to obtain the effective gap equations, in the simple case of the 
FF state. Sec. II I II is devoted to the generalization of this procedure to a generic crystalline 
structure. In Sec. HaZI we show one result of the solution of the smeared gap equations, 
namely the free energy plots for the different crystalline structures. Finally, in Sec. we 
summarize the leading results of our work. 

II. SMEARED LAGRANGIAN IN THE FF PHASE 

We shall consider Cooper pairing of the massless quarks up u and down d, with chemical 
potential fi u , fid- We define /1 = (/x u + //<i)/2 and 5ji = \ji u — yu^|/2 <C ft. We begin with a 
review of the FF state. Although this case can be solved exactly, it is useful to consider it 
here in order to fix the notations and introduce some definitions to be used later on. 

We work in the framework of High Density Effective Theory (HDET) Q Q, H H H 
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29|,|30|. The lagrangian for free quarks can be written as 

C = J2 If+iV ■ dip + + ipUv ■ dtp^ + (L —> R) . (3) 

V 

Here the sum represents an average over velocities; ip± = ip± v are velocity dependent, 
positive energy left handed fields (the negative energy part has been integrated out). ip v 
depends on the residual momentum £, corresponding to the decomposition of the quark 
momentum p — imv + £, with = (0, v) and £» = £-v = £. We also introduce = (1, v) 
and W = (1, -v). 



Next we turn to the interaction term. We consider a Nambu-Jona Lasinio (NJL 
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inspired 



four fermion interaction to mimic the one gluon exchange of QCD, namely 

£ / = _V^ Aa ^7 M A^. (4) 

o 

Here G is a coupling constant, with dimension mass" 2 ; X a are color matrices and a sum over 
flavors is understood (the FF state with quark-quark interaction mediated by one gluon 



exchange has been considered in Ref. 
we get 



In the mean field approximation, after Fierzing, 



d = --e^e^fVf A(r) + c.c.) + (L - R) - -A(r)A*(r) , (5) 

where i, j are flavor indices and a, (3 are color indices. In the FF state the total momentum 
of the Cooper pair is 2q and the condensate has the space-dependence of a single plane 
wave, see Eq. (j2J. In the HDET formalism Eq. (|SJ) can be recast in the form 

£/ = — ^ ^exp{ir- a(vi, vj, q)}e ij e a/ 3 3 ^v i ;i a (x)C^_ Vj;j/3 (x) 

— (L R) + h.c. — -A(r)A*(r) , (6) 
9 

where 

a(vt, vj, q) = 2q - - /ij-vj . (7) 

Eqs. and are the HDET lagrangian of paired quarks in the FF state. Armed with 
them one can write a gap equation, namely a self-consistence (Schwinger-Dyson) equation 
for the gap parameter A. 

Now we have the ingredients necessary to define our smearing procedure. To this end we 
recall the exact FF gap equation [sij ]. 



2jA7iJ 2nJ u £ 2 -e~ K ff 2jAnJ 



where we have defined an effective gap parameter, 



A for (f , v) e PR 
A eff = A9(E u )8(E d ) = { (9) 

elsewhere , 
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33] 



and E U! d are the dispersion laws for u and d quarks respectively 

E d)U = ±6(1 T q • v + + A 2 ; (10) 

PR denotes the pairing region, 

PR= {(f >v ) \E U > O&ndEd > 0} . (11) 

We stress that Eq. (jSJ) is exact. The key observation is that we can obtain the same gap 
equation (jHJ) in the framework of HDET by defining a weighted smearing procedure of the 
gap lagrangian © over the lattice cell. First of all, we note that in the gap equation the 
relevant momenta are small with respect to the gap which is of the order of q. Therefore we 
may assume that the velocity dependent fields are slowly varying over regions of the order 
of the lattice size. This means that in the average we can treat them as constant, and in 
conclusion the average is made only on the coefficient exp{ir ■ a}. Therefore what we are 
computing is 

1(a) = (exp{ir • a}g R {r)} (12) 

where the bracket means average over the cell, and the weight function ^(r) can be chosen 
in such a way that 



/(«) = 6%( |lj , where 5|(x) = j (13) 

elsewhere 

and R/n ~ 1. Independently of the exact form of ^(r), we will assume that the average 
procedure gives as a result the brick-shaped function 5r defined in (113(1 . As shown in 
Ref. 0, 

choosing 

R = , h(v ■ q) = 1 h: , (14) 

2y/e + A 2 |/i(vq)| ' vq' 



for |x| < ^ 
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one obtains from Eq. (JHJ) the smeared lagrangian 

Cj = -\ Y,^effe ij e a03 ^. ia (e)C^ v;jP (-£) - (L - R) + h.c. - ±AA* , (15) 

from which the desired gap eq uation (JSJ) is obtained (the fermion propagator can be read in 
Eqs. (22) and (23) of Ref. jsil). Since R/tt rs 1, then A has not to be small, meaning that 
we should be far from a second order phase transition. 

Let us finally notice that in Eq. (JHJ) one can perform the £q integration by substituting 
the previous expression for R with 

7r|fy-vq| 

i? -2N-Kvq)|- (16) 
In fact, observing that in any case A e yj is equal to or A, at the pole we get back the 
expression (|14|). Then written as in (|16|). i? and analogously A e jj become functions of the 
velocity and the energy; therefore our average should be better taken in the momentum 
space. 



III. SMEARING FOR GENERIC CRYSTALLINE STRUCTURES 

We have defined in the previous section a smearing procedure which, in the case of the 
FF state, allows to write in the HDET formalism a gap equation which coincides with the 
exact gap equation in Eq. (JHJ). Apart the technical details for the definition of the weighting 
function, the smearing procedure can be view as a recipe which allows to replace an in- 
dependent gap function by an r-independent one. All the details of the crystalline structure 
are embodied into the definition of the pairing regions. 

We now use this recipe to smear the interaction lagrangian for generic crystalline struc- 
tures, defined by the pairing ansatz 

p 

A(r) = A ^ e 2iqm r , Qm = q n m . (17) 

m=l 

This procedure, although artificial, allows to study generic crystalline structures in a do- 
main where Ginzburg-Landau expansion may not work well. Generalizing the results of the 
previous equations we substitute in the Lagrangian (JTKJ) A e // (v ■ n, £ ) with 

p 

A E (v, £ ) = A eff (v • n m , £ ) , (18) 

m=l 
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and the arguments of the 9 functions in Eq. Q will be the appropriate fermion dispersion 
laws. After smearing the gap equation reads 

PA-i 9P ! dv ! dH AeM) (19) 
2 J An J 2vr P Q -q-A%(y,to) 1 ' 

which generalizes Eq. (JHJ). The energy integration is performed by the residue theorem and 
the phase space is divided into different regions according to the pole positions. We get 

tff-* , Ae(v,c) tff " Vj - , kA (»> 

where the pairing regions are defined as follows 

ft = {(v,e)|A B (v,e) = *A} (21) 

and we have ruled out the coupling constant G by means of the BCS gap Ao- The first term 
in the sum, corresponding to the region P x , has P equal contributions with a dispersion rule 
equal to the Fulde and Ferrel case. This can be interpreted as a contribution from P non 
interacting plane waves. In the other regions the different plane waves have an overlap. 

Thus for iV plane waves the smearing procedure does not simply leads to a lagrangian 
which is the sum of N plane waves (FF) lagrangians. This "would be nice" scenario is 
complicated by the presence of the pairing regions with two or more overlaps (P2, . . . , Pn)- 



IV. NUMERICAL RESULTS FOR THE FREE ENERGY 

In this section we present results for the free energy Q of the structures that we have 
considered in Ref. jsi|, computed as integral of the gap equation (|2~U|). In particular, in Fig.[T] 
we show fl(A, 5fi) — Q(0, 5 fx) against 5fi for the face centered (fee) and the body centered 
(bec) cubic structures (the other structures have a higher free energy and therefore they are 
not shown here). Our central result is that the bec is a good candidate for the ground state 
for 5fi < 0.95A . Above this value of SfM the fee is the good ground state for 5ji < 1.32A . 
For higher values of 5 /a, the fermion condensation is energetically disfavored. The transition 
to the normal state is first order. 
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FIG. 1: The values of the free energies of the bcc (dashed line) and of the fee (full line) crystalline 
LOFF structures as a function of 5u/Aq. The bcc is the favored structure up to 6fj, ~ .95Ao; for 
.95Ao < Sfj, < 1.32Ao the fee is favored. Here, for each value of 8fi, the values of z q and A are 
those that minimize the free energy. From Ref. 

V. CONCLUSIONS 

In the framework of HDET, we have discussed an averaging procedure of the NJL quark- 
quark interaction lagrangian, treated in the mean field approximation, for the LOFF phase 
of QCD. This procedure gives results which are valid in domains where Ginzburg-Landau 
results may be questionable. Among the several structures considered, we find that a body 
centered cube is the favorite ground state for 5/x < 0.95Ao, while for higher values of 5{i 
the face centered cube is favoured. A first order transition to the normal state is found 
for S/j, ~ 1.32A . The method exposed here can be applied also to the LOFF phase in the 
electromagnetic superconductors (see for example Ref. 
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